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an dem Seil aufgehangten Last (3) mit einem Antriebssystem. Erfindungsgemass weist der Kran oder Bagger eine Bahnsteuerung 
(3 1) auf, deren Ausgangsgrossen direkt oder indirekt als Eingangsgrossen in die Regelung fur die Position bzw. die Geschwindigkeit 
des Krans oder Baggers eingehen, wobei die Fuhrungsgrossen fiirdie Steuerung in der Bahnsteuerung (31) so generiert werden, dass 
sich eine Lastbewegung mit minimierten Pendelausschlagen ergibt. 
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Kran oder Bagger zum Umschlagen von einer an einem Lastseil 
hangenden Last mit optimierter BewegungsfQhrung 



Die Erfindung betrifft einen Kran oder Bagger zum Umschlagen von einer an einem 
Lastseil hangenden Last nach dem Oberbegriff des Anspruchs 1 . 

Im einzelnen. befasst sich die Erfindung mit der Generierung von Fuhrungsgrofien 
als Steuerfunktionen bei Kranen oder Baggern, die eine Bewegung der an einem 
Seil aufgehangten Last in mindestens drei Freiheitsgraden zulasst. Derartige Krane 
oder Bagger weisen ein Drehwerk, das auf einem Fahrwerk aufgebracht sein kann, 
auf, welches zum Drehen des Kranes oder Baggers dient. Weiterhin ist ein Wipp- 
werk zum Aufrichten bzw. Neigen eines Auslegers und ein Schwenkwerk vorhan- 
den. Schlielllich umfasst der Kran oder Bagger ein Hubwerk zum Heben bzw. Sen- 
ken der an dem Seil aufgehangten Last. Derartige Krane oder Bagger finden in ver- 
schiedenster Ausfiihrung Verwendung. Beispielhaft sind hier Hafenmobilkrane, 
Schiffskrane, Offshore-Krane, Raupenkrane bzw. Seilbagger zu nennen. 

Beim Umschlagen einer an einem Seil hangenden Last mittels eines derartigen 
Kranes oder Baggers entstehen Pendelbewegungen der Last, die auf die Bewe- 
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gung des Kranes oder Baggers selbst zurQckzufiihren sind. Es wurden nun bereits 
in der Vergangenheit Anstrengungen unternommen, urn Pendelschwingungen bei 
Lastkranen zu verringern bzw. zu unterdrQcken. 

Die WO 02/32805 A1 beschreibt einen Kran oder Bagger zum Umschlagen von 
einer an einem Lastseil hangenden Last mit einer computergesteuerten Regelung 
zur Dampfung der Lastpendelung, die ein Bahnplanungsmodul, eine Zentripetal- 
kraftkompensationseinrichtung und zumindest einen Achsregler fur das Drehwerk, 
einen Achsregler fur das Wippwerk und einen Achsregler fur das Hubwerk aufweist. 
Dabei wird im Bahnplanungsmodul nur die kinematischen Begrenzungen des Sys- 
tems beriicksichtigt. Das dynamische Verhalten findet nur beim Entwurf der Rege- 
lung Berucksichtigung. 

Aufgabe der Erfindung ist es, die BewegungsfQhrung der am Lastseil hangenden 
Last noch weiter zu opimieren. 

Zur LSsung dieser Aufgabe weist ein gattungsgemalier Kran oder Bagger eine 
Steuerung auf, in der die FuhrungsgroRen fQr die Steuerung so generiert werden, 
dass sich eine optimierte Bewegung mit minimierten Pendelausschlagen ergibt. 
Hierbei kann auch die abgefahrene Bahn der pendelnden Last prognostiziert wer- 
den und darauf aufbauend eine Kollisionsvermeidungsstrategie realisiert werden. 

Vorteilhafte Ausgestaltungen der Erfindung ergeben sich aus den sich an den 
Hauptanspruch anschliedenden UnteransprQchen. 

Von besonderem Vorteil ist es, dass in der Bahnsteuerung der vorliegenden Erfin- 
dung modelbasiert optimale Steuertrajektorien online berechnet und aktualisiert 
werden. Dabei k6nnen die modelbasierten optimalen Steuertrajektorien basierend 
auf einem urn Referenztrajektorien linearisierten Model erstellbar sein. Alternativ 
kSnnen die modelbasierten optimalen Steuertrajektorien auf einem nicht linearen 
Modelansatz basieren. 
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Die modelbasierten optimalen Steuertrjektorien konnen unter ROckfuhrung von al- 
ien Zustandsgroften ermittelt werden. 

Alternativ konnen die modelbasierten optimalen Steuertrajektorien unter RQckfuh- 
rung von mindestens einer Messgrofce und Schatzung der verbleibenden Zu- 
standsgrdfien ermittelbar sein. 

Wiederum alternativ kdnnen die modelbasierten optimalen Steuertrajektorien unter 
Ruckfuhrung von mindestens einer MessgrClie und Nachfuhrung der verbleibenden 
Zustandsgrofcen durch modelbasierte Vorwartssteuerung ermittelt werden. 

Die Bahnsteuerung kann vorteilhaft als Vollautomatik oder aber auch als Halbau- 
tomatik durchfuhrbar sein. 

So ergibt sich im Zusammenhang mit einer Regelung zur Lastpendeldampfung ein 
optimiertes Bewegungsverhalten mit vermindertem Restpendeln und geringeren 
Pendelausschlagen wahrend der Fahrt. Ohne die Regelung zur Lastpendeldamp- 
fung kann die erforderliche Sensorik am Kran verringert werden. Es kann ein voll- 
automatischer Betrieb, bei dem Start und Zielpunkt feststehen ebenso realisiert 
werden, wie ein Handhebelbetrieb, der im folgenden als halbautomatischer Betrieb 
bezeichnet wird. 

In der vorliegenden Erfindung werden die Sollfunktionen im Gegensatz zur WO 
02/32805 A1 nun derart generiert, dass bereits vor dem Aufschalten auf die Rege- 
lung das dynamische Verhalten des Krans berucksichtigt ist. Damit hat die Rege- 
lung nur noch die Aufgabe, Modellabweichungen und StOrgroBen auszugleichen, 
wodurch ein verbessertes Fahrverhalten resultiert. Zudem kann, wenn die Positi- 
onsgenauigkeit und das tolerierbare Restpendeln es zulasst, die Regelung ganz 
entfallen und der Kran mit dieser optimierten Steuerfunktion betrieben werden. Al- 
lerdings wird das Verhalten etwas ungunstiger als beim Betrieb mit der Regelung 
sein, da das Modell nicht in alien Details mit den tatsachlichen Gegebenheiten Q- 
bereinstimmt. 
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Das Verfahren sieht zwei Betriebsmodi vor. Den Handhebelbetrieb, bei dem der 
Bediener durch die Handhebelauslenkung eine Sollgeschwindigkeit der Last vor- 
gibt, und dem vollautomatischen Betrieb, bei dem Start- und Zielpunkt vorgegeben 
werden. 

Daneben kann die optimierte Steuerfunktionsberechnung allein oder in Zusammen- 
hang mit einer Regelung zur Lastpendeldampfung betrieben werden. 

Weitere Einzelheiten und Vorteile der Erfindung werden anhand eines in der Zeich- 
nung dargestellten AusfQhrungsbeispiels erlautert. Als typischer Vertreter fur einen 
Kran oder Bagger der eingangs genannten Gattung wird die Erfindung hier anhand 
eines Hafenmobilkranes beschrieben. 

Weitere Einzelheiten und Vorteile der Erfindung werden anhand eines in der Zeich- 
nung dargestellten AusfQhrungsbeispiels erlautert. Als typischer Vertreter fur einen 
Kran oder Bagger der eingangs genannten Gattung wird die Erfindung hier anhand 
eines Hafenmobilkranes beschrieben. 

Es zeigen: 

Fig. 1 : Prinzipielle mechanische Struktur eines Hafenmobilkranes 

Fig. 2: Zusammenwirken von hydraulischer Steuerung und Bahnsteuerung mit 
Modul zur optimierten BewegungsfUhrung als Steuerfunktion des Kranes 

Fig. 3: Struktur der Bahnsteuerung mit Modul zur optimierten BewegungsfOh- 
rung mit Regelung zur Lastpendeldampfung 

Fig. 4: Struktur der Bahnsteuerung mit Modul zur optimierten Bewegungsfuh- 
rung als Steuerfunktion ohne Regelung zur Lastpendeldampfung (ggf. 
mit unterlagerten Positionsreglern fur die Antriebe) 

Fig. 5: Mechanischer Aufbau des Drehwerks und Definition von Modellvariablen 

Fig. 6: Mechanischer Aufbau des Wippwerks und Definition von Modellvariablen 

Fig. 7: Aufrichtkinematik des Wippwerks 
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Fig. 8: Ablaufdiagramm fur die Berechnung der optimierten Steuergrofie im 

vollautomatischen Betrieb 
Fig. 9: Ablaufdiagramm fQr die Berechnung der optimierten Steuergrolie im 

halbautomatischen Betrieb 
Fig. 10: Beispielhafte FuhrungsgrGBengenerierung im vollautomatischen Betrieb 
Fig. 1 1 : Beispielhafte Zeitverlaufe von SteuergrSlien und Regelgrofien im Hand- 

hebelbetrieb 

In Fig. 1. ist die prinzipielle mechanische Struktur eines Hafenmobilkrans darge- 
stellt. Der Hafenmobilkran ist zumeist auf einem Fahrgestell 1 montiert. Zur Positio- 
nierung der Last 3 im Arbeitsraum kann der Ausleger 5 mit dem Hydraulikzylinder 
des Wippwerks 7 urn den Winkel <p A gekippt werden. Mit dem Hubwerk kann die 
Seillange / s variiert werden. Der Turm 1 1 ermoglicht die Drehung des Auslegers um 
den Winkel <po um die Hochachse. Mit dem Lastschwenkwerk 9 kann die Last am 
Zielpunkt um den Winkel <p mt gedreht werden. 

Fig. 2 zeigt das Zusammenwirken von hydraulischer Steuerung und Bahnsteuerung 
31 mit Modul zur optimierten Bewegungfuhrung. In der Regel besitzt der Hafenmo- 
bilkran ein hydraulisches Antriebssystem 21. Ein Verbrennungsmotor 23 speist 0- 
ber ein Verteilergetriebe die hydraulischen Steuerkreise. Die hydraulischen Steuer- 
kreise bestehen jeweils aus einer Verstellpumpe 25, die uber ein Proportionalventil 
im Vorsteuerkreis angesteuert wird, und einem Motor 27 oder Zylinder 29 als Ar- 
beitsmaschine. Ober das Proportionalventil wird damit lastdruckunabhangig ein 
Forderstrom Qfd, Qfa, Qfu Qfr eingestellt. Die Proportionalventile werden Qber die 
Signale u S to, ustA, u S tu ustR angesteuert. Die hydraulische Steuerung ist meist mit 
einer unterlagerten Forderstromregelung ausgestattet. Wesentlich ist dabei, dalJ die 
Steuerspannungen u S tD, "sm, Ustu u^r an den Proportionalventilen durch die un- 
terlagerte Forderstromregelung in hierzu proportionale Forderstrome Qfd, Qfa, 
Qfl, Qfr im entsprechenden Hydraulikkreislauf umgesetzt werden. 

Die Struktur der Bahnregelung ist nun in den Figuren 3 und 4 dargestellt. Figur 3 
zeigt die Bahnsteuerung mit dem Modul zur optimierten Bewegungsfiihrung mit 
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Regelung zur Lastpendeldampfung und Figur 4 die Bahnsteuerung mit dem Modul 
zur optimierten Bewegungsfuhrung ohne Regelung zur Lastpendeldampfung. Diese 
Lastpendeldampfung kann beispielsweise nach der Schrift PCT/EP01/12080 ent- 
worfen worden sein. Daher wird der dort offenbarte Inhalt vollumfanglich in diese 
Schrift mit einbezogen. 

Wesentlich ist nun, daft die Zeitfunktionen fur die Steuerspannungen der Proportio- 
nalventile nicht mehr direkt aus den Handhebeln beispielsweise Qber Rampenfunk- 
tionen oder einen Bahnplaner, der die kinematischen Beschrankungen des System 
beriicksichtigt, abgeleitet werden, sondern derart in der Bahnsteuerung 31 berech- 
net werden, daft beim Bewegen des Krans keine oder geringe Pendelbewegungen 
der Last auftreten und die Last der gewQnschten Bahn im Arbeitsraum folgt. D.h. 
bei der Berechnung der optimierten Steuergrolie wird nicht nur die kinematische 
Beschreibung sondern die dynamische Beschreibung des Systems beriicksichtigt. 

EingangsgrSfien des Moduls 37 ist eine Sollpunktmatrix 35 fur die Position und Ori- 
entierung der Last, die im einfachsten Fall aus Start- und Zielpunkt besteht,. Die 
Position wird ublicherweise bei Drehkranen durch Polarkoordinaten beschrieben 
(<Pld, riA, I). Da dadurch die Position des ausgedehneten KSrpers (bspw. eines 
Containers) im Raum nicht vollstandig beschrieben wird, kann eine weitere Winkel- 
grofte hinzugenommen werden (Drehwinkel y L urn die Hochachse, die parallel zum 
Seil liegt). Die ZielpositionsgrOflen <p LD ztei, runei, Izia , nztei sind im Vektor qziei zu- 
sammengefasst. 

Eingangsgroften des Moduls 39 sind die aktuellen Stellungen der Handhebel 34 zur 
Ansteuerung des Krans. Die Auslenkung der Handhebel entspricht der gewQnsch- 
ten Zielgeschwindigkeit der Last in der jeweiligen Bewegungsrichtung. Dement- 
sprechend werden die Ziel-Geschwindigkeiten <$j J DZieU i 'LAZieijzieiAziei' ZX3X x\ Ziel- 
geschwindigkeitsvektor q zusammengefasst. 



WO 2004/106215 



-7- 



PCT/EP2004/005734 



lm Falle des Moduls zur optimierten BewegungsfOhrung im vollautomatischen Be- 
trieb 37 kann aus dieser Information uber das abgelegte Modell zur Beschreibung 
des dynamischen Verhaltens und den gewahlten Rand- und Nebenbedingungen 
das Optimalsteuerungsproblem gelost werden. Ausgangsgrofien sind dann die 
Zeitfunktionen u ollttD , u out ,A> u outt u u outtR , die zugleich EingangsgrSften der unterlager- 
ten Regelung zur Lastpendeldampfung 36 bzw. der unterlagerten Regelung fur Po- 
sition bzw. Geschwindigkeit des Krans 41. Auch eine direkte Ansteuerung 41 des 
Krans ohne unterlagerte Regelung ist bei entsprechender Formulierung der Glei- 
chungen in 37 moglich. Dabei kann im vollautomatischen Betrieb der Handhebel- 
wert dazu benutzt werden, im Optimalsteuerungsproblem die Nebendbedingung 
der maximalen zulassigen Geschiwndigkeit zu verandern. Dies ist insbesondere 
deshalb vorteilhaft, daft auch im vollautomatischen Betrieb der Anwender die Mog- 
lichkeit hat, den vollautomatischen Ablauf online in der Geschwindigkeit zu beein- 
flussen. Die vorgenommenen Anderungen werden sofort im nachsten Durchlauf 
des Algorithmus ubernommen und berucksichtigt. 

Im Falle des Moduls zur optimierten Beweungsfuhrung im halbautomatischen Be- 
trieb 39 wird jedoch zur Information der aktuell gewunschten Zielgeschwindigkeit 
der Last durch die Handhebelstellung als weitere Information der aktuelle System- 
zustand neben der Rand- und Nebenbedingungen benotigt. Deshalb miissen im 
halbautomatischen Betrieb die Meligroften der Position von Kran und Last laufend 
auf das Modul 39 zuruckgefQhrt werden. Im einzelnen sind dies: 

- Drehwerkswinkel <p D , 

- Wippwerkswinkel <p A , 

- Seillange l s , und 

- relative Lasthakenposition c 

und die Winkel zur Beschreibung der Lastposition: 

- tangentialer Seilwinkel 

- radialer Seilwinkel <p Sn und 



WO 2004/106215 



PCT/EP2004/005734 



-8- 

- absoluter Rotationswinkel der Last y L . 

Insbesondere letztgenannte MefSgroften fur Seilwinkel und absoluten Rotationswin- 
kel der Last sind nur mit grolierem Aufwand messtechnisch zu erfassen. Fur die 
Realisierung einer Lastpendeldampfung sind diese jedoch unabdingbar notwendig, 
urn Storungen auszugieichen. Dadurch kann eine sehr hohe Positioniergenauigkeit 
bei geringem Restpendeln auch unter EinfluB von StorgrOBen (wie Wind) erreicht 
werden. Im Falle von Fig. 3 stehen diese Gr6Ben alle zur VerfQgung. 

Wird jedoch das Verfahren in einem System eingesetzt, in dem keine Sensoren fQr 
die Seilwinkelmessung und den absoluten Rotationswinkel existieren, so muB fur 
das Modul zur optimierten Bewegungsfuhrung im halbautomatischen Betrieb diese 
GroBen rekonstruiert werden. Hier bieten sich modellbasierte Schatzverfahren 43, 
wie Beobachterstrukturen, an. Hierbei wird aus den MeBgroBen der Kranposition 
und den Ansteuerfunktionen u otlltD , u ou ,,a, «<>«/,/> in einem hinterlegten dynami- 
schen Modell die fehlenden Zustandsgrofcen geschatzt oder nachgefuhrt (siehe Fig. 
4). 

Grundlage fur das Verfahren der optimierten Bewegungsfuhrung ist das Verfahren 
der dynamischen Optimierung. Hierzu mu(i das dynamische Verhalten des Kranes 
in einem Differentialgleichungsmodell abgebildet werden. Zur Ableitung der Modell- 
gleichungen kann entweder der Lagrange-Formalismus oder die Methode nach 
Newton-Euler verwendet werden. 

Im folgenden werden mehrere mogliche Modellansatze vorgestellt. Zunachst wer- 
den anhand Fig. 5 und 6 die Definition der Modellvariablen vorgenommen. Zur bes- 
seren Obersichtlichkeit zeigt Fig. 5 die Modellvariablen die im Zusammenhang mit 
der Drehbewegung stehenden Modellvariablen und Fig. 6 die Modellvariablen fQr 
die radiale Bewegung. 

Zunachst wird Fig. 5 detailliert erlautert. Wesentlich ist dabei der dort gezeigte Zu- 
sammenhang zwischen der Drehposition q>o des Kranturmes und der Lastposition 
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<p LD in Drehrichtung. Die urn den Pendelwinkel korrigierte Lastdrehwinkelposition 
berechnet sich dann zu 



Is ist dabei die resultierende SeillSnge vom Auslegerkopf bis zum Lastmittelpunkt. 
<p A ist der aktuelle Aufrichtwinkel des Wippwerks, l A ist die Lange des Auslegers, <p S t 
ist der aktuelle Seilwinkel in tangentialer Richtung (da tpst klein ist, kann sm<p S t*ipst 
angenahert werden). 

Das dynamische System fur die Bewegung der Last in Drehrichtung kann durch die 
folgenden Differentialgleichungen beschrieben werden. 

[j T +(JAZ+ m A s A +m L l%)cos 2 <p A \p D +m L l A l s cosq> A cp st +b D <p D =M MD -Af^(2) 
m L l A l s COS( PA 9D + m L Is <Pst + m LS h <Pst = 0 ( 3 ) 

Bezeichnungen: 

m L Lastmasse 

l s Seillange 

m A Masse des Auslegers 

Jaz Massentragheitsmoment des Auslegers bezuglich Schwerpunkt bei 



<Pu> =< Pd+ arctan 



l s9st 



d) 



I A COS( PA 



Jt 
b D 
Mmd 
M rd 



Drehung um Hochachse 
Lange des Auslegers 
Schwerpunktsabstand des Auslegers 
Massentragheitsmoment des Turmes 
viskose Dampfung im Antrieb 
Antriebsmoment 



Reibmoment 
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(2) beschreibt im wesentlichen die Bewegungsgleichung fQr den Kranturm mit Aus- 
leger, wobei die RQckwirkung durch die Lastpendelung berficksichtigt wird. (3) ist 
die Bewegungsgleichung, welche die Lastpendelung urn den Winkel <pst beschreibt, 
wobei die Anregung der Lastpendelung durch die Drehung des Turmes Qber die 
Winkelbeschleunigung des Turmes oder eine aufSere Stdrung, ausgedrtickt durch 
Anfangsbedingungen fOr diese Differentialgleichungen, verursacht wird. 

Der hydraulische Antrieb wird durch die folgenden Gleichungen beschrieben. 



V 

M M D=iD^ A PD 

1 V 

**D = ^(QfD ~ iD ) ( 4 > 

QfD ~ K PD u StD 



i D ist das Obersetzungsverhaltnis zwischen Motordrehzahl und Drehgeschwindig- 
keit des Turms, V ist das Schluckvolumen der Hydraulikmotoren, Apo ist der Druck- 
abfall Qber dem hydraulischen Antriebsmotor, /? ist die Olkompressibilitat, Qfd ist 
der Forderstrom im Hydraulikkreis fQr das Drehen und K PD ist die Proportionalitats- 
konstante, die den Zusammenhang zwischen Forderstrom und Ansteuerspannung 
des Proportionalventils angibt. Dynamische Effekte der unterlagerten Forderstrom- 
regelung werden vernachlassigt. 

Alterriativ hierzu kann das Obertragungsverhalten der Antriebsaggregate anstatt mit 
der Gleichung 4 durch einen approximativen Zusammenhang als Verzogerungs- 
glied 1. oder hSherer Ordnung dargestellt werden. Im folgenden ist die Approxima- 
tion mit einem Verzogerungsglied 1. Ordnung dargestellt. Danach ergibt sich die 
Obertragungsfunktion 



s® D (s)= U slD (s) (5) 



bzw. im Zeitbereich 
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T, 



DAntr 



1 




U 



StD 



(6) 



Damit kann aus den Gleichungen (6) und (3) ebenfalls eine adaquate Modellbe- 

schreibung aufgebaut werden; Gleichung (2) wird nicht ben6tigt. 

T D Antr ist die approximative (aus Messungen ermittelte Zeitkonstante zur Beschrei- 

bung des Verzogerungsverhaltens der Antriebe. K PDAnt r die resultiertende Verstar- 

kung zwischen Ansteuerspannung un resultierender Geschwindigkeit im staionaren 

Fall. 

Bei einer unerheblichen Zeitkonstante bezuglich der Antriebsdynamik kann direkt 
eine Proportionalitat zwischen Geschwindigkeit und Ansteuerspannung des Pro- 
portionalventils angenommen werden. 

<PD= K PDdirekt u StD i 7 ) . 

Auch hier kann dann aus den Gleichungen (7) und (3) eine adaquate Modellbe- 
schreibung aufgebaut werden. 

Fur die in Fig. 6 dargestellte radiale Bewegung lasst sich analog zu den Gleichun- 
gen (2) und (3) die Bewegungsgleichungen aufstellen. Hierzu gibt Fig. 6 Erlaute- 
rungen zur Definition der Modellvariablen. Wesentlich ist dabei der dort gezeigte 
Zusammenhang zwischen der Aufrichtwinkelposition <p A des Auslegers und der 
Lastposition in radialer Richtung rtA 



r LA =Ia COS( Pa +l S<Psr 



(8) 



Das dynamische System kann dann nach Anwendung des Newton-Euler- 
verfahrens durch die folgenden Differentialgleichungen beschrieben werden. 
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[j AY + m A s A + m L l\ sin 2 <p A Jis) A - m L l A l s sin9^ ip sr 

+ b A q> A - m A s A g sin<p^ • 9^ = (9) 
M MA -Mra -m A s A g coses? A 

~ m L l A l s sin( P^ <?A + m L l s Vsr + m L l s S 9sr = m L l s $D Vs9sr +l A C( > s 9^) 

(10) 

Bezeichnungen: 



nriL Lastmasse 

/ s Seillange 

itia Masse des Auslegers 

J AY Massentragheitsmoment bezQglich Schwerpunkt bei Drehung urn 

horizontale Achse inkl. Antriebsstrang 

l A Lange des Auslegers 

sa Schwerpunktsabstand des Auslegers 

bA viskose Dampfung 

Mma Antriebsmoment 



Mra Reibmoment 

Gleichung (9) beschreibt im wesentlichen die Bewegungsgleichung des Auslegers 
mit dem antreibenden Hydraulikzylinder, wobei die Ruckwirkung durch die Pende- 
lung der Last beriicksichtigt wird. Dabei ist auch der durch die Schwerkraft des 
Auslegers einwirkende Anteil und die viskose Reibung im Antrieb berOcksichtigt. 
Gleichung (10) ist die Bewegungsgleichung, welche die Lastpendelung <p Sr be- 
schreibt, wobei die Anregung der Schwingung durch das Aufrichten bzw. Neigen 
des Auslegers Ober die Winkelbeschleunigung des Auslegers Oder eine aufcere 
Storung, ausgedrOckt durch Anfangsbedingungen for diese Differentialgleichungen, 
verursacht wird. Ober den Term auf der rechten Seite der Differentialgleichung wird 
der EinfluB der Zentripetalkraft auf die Last bei Drehung der Last mit dem Drehwerk 
beschrieben. Dadurch wird ein ftir einen Drehkran typisches Problem beschrieben, 
da damit eine Kopplung zwischen Drehwerk und Wippwerk besteht. Anschaulich 
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kann man dieses Problem dadurch beschreiben, daft eine Drehwerksbewegung mit 
quadratischer Drehgeschwindigkeitsabhangigkeit auch einen Winkelausschlag in 
radialer Richtung hervorruft. 

Der hydraulische Antrieb wird durch die folgenden Gleichungen beschrieben. 



M MA = F Zyl d b cos <p p (<p A ) 
F Zyl = PZyl A Zyl 

2 (11) 

PZyl = — (QFA ~ A Zyl z Zyl (<PA ><PA )) 
P v Zyl 

QFA = K PA u StA 



F Zy i ist die Kraft des Hydraulikzylinders auf die Kolbenstange, p Zy i ist der Druck im 
Zylinder (je nach Bewegungsrichtung kolben- oder ringseitig), A Zy i ist die Quer- 
schnittsflache des Zylinders Qe nach Bewegungsrichtung kolben- oder ringseitig), ft 
ist die Olkompressibilitat, V Zy i ist das Zylindervolumen, Q FA ist der Forderstrom im 
Hydraulikkreis fur das Wippwerk und K PA ist die Proportionalitatskonstante, die den 
Zusammenhang zwischen Forderstrom und Ansteuerspannung des Proportional- 
ventils angibt. Dynamische Effekte der unterlagerten Forderstromregelung werden 
vemachlSssigt. Bei der Olkompression im Zylinder wird als relevantes Zylindervo- 
lumen die Halfte des Gesamtvolumens des Hydraulikzylinders angenommen. 
z Zy i,z Zy i si nd die Position bzw. die Geschwindigkeit der Zylinderstange. Diese sind 

ebenso wie die geometrischen Parameter d/> und <p p von der Aufrichtkinematik ab- 
hangig. 

In Fig. 7 ist die Aufrichtkinematik des Wippwerks dargestellt. Beispielhaft ist der 
Hydraulikzylinder oberhalb des Drehpunktes des Auslegers am Kranturms veran- 
kert. Aus Konstruktionsdaten kann der Abstand d a zwischen diesem Punkt und dem 
Drehpunkt des Auslegers entnommen werden. Die Kolbenstange des Hydraulikzy- 
linders ist am Ausleger im Abstand cfe befestigt. Der Korrekturwinkel <po berOcksich- 
tigt die Abweichungen der Befestigungspunkte von der Ausleger- bzw. Turmachse 
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und ist ebenfalls aus Konstruktionsdaten bekannt. Daraus lafct sich der folgende 
Zusammenhang zwischen Aufrichtwinkel <p A und Hydraulikzylinderposition z Zy / her- 
leiten. 

z Zyl = 4d] + d\ - 2d b d a sm(<p A - <p 0 ) (1 2) 

Da nur der Aufrichtwinkel <p A Mefcgrofce ist, ist die umgekehrte Relation von (12) 
sowie die Abhangigkeit zwischen Kolbenstangengeschwindigkeit z Zyl und Aufricht- 

geschwindigkeit ^ebenfalls von Interesse. 



da d b — z\ t 

<p A = arcsin( —— ^) + <p Q (1 3) 

2d a d b 



. _ dq> A . _ 4 d l + d] - 2d b d a sin(^ - <p Q ) . 



dz Zyl 9 -d b d a cos(<p A -<p 0 ) 



Fur die Berechnung des wirksamen Momentes auf den Ausleger ist aufierdem die 
Berechnung des Projektionswlnkels <p p erforderlich. 

cosp,= . (15) 
V rf « + rf * ~ 2rf W* sin(^ - <p 0 ) 

Altemativ kann hierzu wieder anstatt der Hydraulikgleichungen (11) eine NSherung 
fur die Dynamik der Antriebe mit einen approximativen Zusammenhang als Verzo- 
gerungsglied 1. oder hoherer Ordnung vorgesehen werden. Damit erhalt man bei- 
spielhaft 

K c 



^w s u r u M t*) (16) 



^PAAntr 
^ AAntr 

bzw. im Zeitbereich 



S*--^***^**-"** (17) 

* AAntr 1 AAntr 
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Damit kann aus den Gleichungen (17), (14) und (10) ebenfalls eine adaquate Mo- 
dellbeschreibung aufgebaut werden; Gleichung (9) wird nicht benotigt. T AAn tr ist die 
approximative (aus Messungen ermittelte Zeitkonstante zur Beschreibung des Ver- 
zSgerungsverhaltens der Antriebe. K PA Antr die resultiertende Verstarkung zwischen 
Ansteuerspannung und resultierender Geschwindigkeit im stationaren Fall. 
Bei einer unerheblichen Zeitkonstante bezuglich der Antriebsdynamik kann direkt 
eine Proportionalitat zwischen Geschwindigkeit und Ansteuerspannung des Pro- 
portionalventils angenommen werden. 

ZZyl = K PAdirekt u StA ( 1 8 ) 

Auch hier kann dann aus den Gleichungen (18), (10) und (14) eine adaquate Mo- 
dellbeschreibung aufgebaut werden. 

Letzte Bewegungsrichtung ist das Drehen der Last am Lasthaken selbst durch das 
Lastschwenkwerk. Eine entsprechende Beschreibung dieser Regelung ergibt sich 
aus der deutschen Patentanmeldung DE 100 29 579 vom 15.06.2000, auf deren 
Inhalt hier ausdrOcklich verwiesen wird. Die Rotation der Last wird Uber das zwi- 
schen einer am Seil hangenden Unterflasche und einer Lastaufnahmevorrichtung 
angeordnete Lastschwenkwerk vorgenommen. Dabei werden auftretende Torsi- 
onsschwingungen unterdruckt. Damit kann die in den meisten Fallen ja eben nicht 
rotationssymmetrische Last lagegenau aufgenommen, durch einen entsprechenden 
Engpass bewegt und abgesetzt werden. Selbstverstandlich wird auch diese Bewe- 
gungsrichtung im Modul zur optimierten BewegungsfOhrung integriert, wie dies bei- 
spielsweise anhand der Obersicht in Fig. 3 dargestellt ist. In besonders vorteilhafter 
Weise kann hier die Last schon nach dem Aufnehmen wahrend des Transportes 
durch die Luft in die entsprechend gewunschte Schwenkposition mittels des Last- 
schwenkwerkes verfahren werden, wobei hier die einzelnen Pumpen und Motoren 
synchron angesteuert werden. Wahlweise kann auch ein Modus fOr eine drehwin- 
kelunabhangige Orientierung gewahlt werden. 
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Daraus ergibt sich die nachfolgend aufgefQhrte Bewegungsgleichung. Die Varai- 
ablenbezeichnung entsprechen der DE 100 29 579 vom 15.06.2000. Dabei wurde 
keine Linearisierung vorgenommen. 

(0^ +® Uc )y d mi =-m L gsin(^^-)^--e Lc c (19) 

Auch fOr das Lastschwenkwerk konnen nun Differentialgleichungen zur Beschrei- 
bung der Antriebsdynamik zur Verbesserung der Funktion wie bei der Drehbewe- 
gung zusatzlich berucksichtigt werden. Hier soil auf eine detaillierte Darstellung 
verzichtet werden. 

Die Dynamik des Hubwerks sei vernachlassigt, da die Dynamik der Hubwerksbe- 
wegung im Vergleich zur Systemdynamik des Lastpendelung des Krans schnell ist. 
Wie beim Lastschwenkwerk konnen jedoch bei Bedarf die entsprechenden dynami- 
schen Gleichungen zur Beschreibung der Hubwerksdynamik jederzeit erganzt wer- 
den. 

Die verbleibenden Gleichungen zur Beschreibung des Systemverhaltens sollen nun 
in eine nichtlineare Zustandsraumdarstellung nach Isidori, Nonlinear Control Sys- 
tems Springer Verlag 1995 gebracht werden. Dies sei beispielhaft basierend auf 
den Gleichungen (2), (3), (9), (10), (14), (15) durchgefOhrt. Dabei ist die Rotations- 
achse der Last urn die Hochachse und die Hubwerksachse in diesem nun nachfol- 
genden Beispiel nicht berucksichtigt. Es ist jedoch keine Schwierigkeit, diese in die 
Modellbeschreibung mit aufzunehmen. FQr den vortiegenden Andwendungsfall sei 
ein Kran ohne automatisches Lastschwenkwerk angenommen, das Hubwerk wird 
aus SicherheitsgrOnden vom Kranbediener manuell bedient. Dementsprechend er- 
halt man: 

x = a(x) + k(x)u 

Zustandsraumdarstellung: . . (20) 

y = c(x) 



mit: 
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Zustandsvektor: 



^ = [9^ <p D q>A VA <?St <PSr <PSr <PSr PZyl 



(21) 



SteuergroBe: 



u = [u StD u SfA ] T 



(22) 



AusgangsgroBe: 



y = b>LD r LA \ 



(23) 



Die Vektoren a(x),6(x),c(x)ergeben sich durch Umformung der Gleichungen (2)- 
(4), (8)-(15). 

Im Betrieb des Moduls zur optimierten Bewegungsfuhrung ohne unterlagerte Last- 
pendeldampfung tritt im halbautomatischen Betrieb die Problematik auf, dass der 
Zustand x als Messvektor vollstandig vorliegen muB. Da in diesem Fall aber keine 
Pendelwinkelsensoren installiert sind, mQssen die in diesem oben beschriebenen 
Fall beispielhaft die PendelwinkelgroBen <vst>yst>WSrAsr aus den AnsteuergroBen 
ustD, u s ,a und den MeBgroBen <s> D ,<?-D><?A><f>A>Pzyl rekonstruiert werden. Hierzu 
wird das nichtlineare Modell nach Gleichung (20-23) linearisiert und beispielsweise 
ein parameteradaptiver Zustandsbeobachter (siehe auch Fig. 4 Block 43) entwor- 
fen. Bei reduzierten Genauigkeitsanforderungen kann vereinfachend auch eine Zu- 
standsnachfuhrung der SeilwinkelgrSBen basierend auf den Modellgleichungen und 
den bekannten Verlaufen der EingangsgroBen sowie der messbaren Zustandsgro- 
Ben erfolgen. 

Die Sollverlaufe fOr die Eingangssignaie (SteuergrSBen) u s<o (0,«^(0 werden 
durch die Losung eines Optimalsteuerungsproblems, d.Dh. einer Aufgabe der dy- 
namischen Optimierung bestimmt. Hierzu wird die angestrebte Reduktion des Last- 
pendelns in einem Zielfunktional erfasst. Randbedingungen und Trajektorien- 
beschrankungen des Optimalsteuerungsproblems ergeben sich aus den Bahnda- 
ten, den technischen Restriktionen des Kransystems (z.DB. limitierte Antriebsleis- 
tung, sowie Beschrankungen aufgrund von dynamischer Lastmomentbegrenzungen 
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zur Verhinderung des Kippens des Krans) sowte erweiterten Forderungen an die 
Bewegung der Last. Beispielsweise ist es mit dem im folgenden beschriebenen 
Verfahren erstmals moglich, den Bahnkorridor, den die Last bei der Aufschaltung 
der berechneten Steuerfunktionen bendtigt, exakt im vorhinein zu pradizieren. Da- 
mit sind Automatisierungsmoglichkeiten gegeben, die vormals nicht losbar waren. 
Eine solche Formulierung des Optimalsteuerungsproblems wird im folgenden bei- 
spielhaft sowohl fQr den vollautomatischen Betrieb des Systems mit vorgegebenem 
Start- und Zielpunkt der Lastbahn als auch fur den Handhebelbetrieb gegeben. 

Im Fall des vollautomatischen Betriebs wird die gesamte Bewegung vom vorgege- 
benem Start- bis zum vorgegebenen Zielpunkt betrachtet. Im Zielfunktional des Op- 
timalsteuerungsproblems werden die Lastpendelwinkel quadratisch bewertet. Die 
Minimierung dieses Zielfunktionals liefert daher eine Bewegung mit reduzierter 
Lastpendelung. Eine zusatzliche Bewertung der Lastpendelwinkelgeschwindigkei- 
ten mit einem zeitvarianten (zum Ende des Optimierungshorizonts zunehmenden) 
Strafterm ergibt eine Beruhigung der Lastbewegung am Ende des Optimierungsho- 
rizonts. Ein Regularisierungsterm mit quadratischer Bewertung der Amplituden der 
SteuergroBen kann die numerische Kondition der Aufgabe gunstig beeinflussen. 

J = )(<Pl (0 + <Pl (f) + P(0fe (0 + 9l (0)+ Pu (Ms,d (0» "sm (0))dt (24) 

'» 

Bezeichnungen: 

t 0 Vorgegebener Startzeitpunkt 

t f Vorgegebener Endzeitpunkt 

p(t) Zeitvarianter Strafkoeffizient 

Pu(ust D ,UstA) Regularisierungsterm (quadratische Bewertung der Steuergrblien) 

Im Handhebelbetrieb wird dagegen nicht die komplette Lastbewegung zwischen 
vorgegebenem Start- und Zielpunkt betrachtet, sondem das Optimalsteuerungs- 
problem wird auf einem mit dem dynamischen Vorgang mitbewegten Zeitfenster 
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[t 0 J f ] betrachtet. Der Startzeitpunkt des Optimierungshorizonts * 0 ist der aktuelle 

Zeitpunkt, und im Optimalsteuerungsproblem wird die Dynamik des Kransystems 
im Prognosehorizont bis t f betrachtet. Dieser Zeithorizont ist ein wesentlicher Tu- 
ning-Parameter des Verfahrens und ist nach unten durch die Schwingungsdauer 
der Lastpendelbewegung begrenzt. 

Im Zielfunktional des Optimalsteuerungsproblems ist neben der angestrebten Re- 
duktion des Lastpendelns die Abweichung der tatsachlichen Lastgeschwindigkeit 
von den durch die Handhebelsteilungen vorgegebenen Sollgeschwindigkeiten zu 
berQcksichtigen. 

f\+ VlQ) + <Plr(!) + P(0fe (0 + *l(0)+ Pu (««D <f\u su (0)J ^ ° } 

t 

Bezeichnungen: 



t 0 Vorgegebener Startzeitpunkt des Optimierungshorizonts 

t f Vorgegebener Endzeitpunkt des Prognosezeitraums 

Pld Bewertungskoeffizient Abweichung Lastdrehwinkelgeschwindigkeit 

Pld** Durch Handhebelstellung vorgegebene Lastdrehwinkelgeschwin- 

digkeit 

Pla Bewertungskoeffizient Abweichung radiale Lastgeschwindigkeit 

i'LAsou Durch Handhebelstellung vorgegebene radiale Lastgeschwindig- 

keit 



Im vollautomatischen Betrieb mit vorgegebenem Start- und Zielpunkt ergeben sich 
die Randbedingungen far das Optimalsteuerungsproblem aus deren Koordinaten 
und den Anforderungen einer Ruhelage in Start und Zielposition. 
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<Pd Co ) = ^d.o » <Pt> (ff ) = <Pdj 

<j>M = o, <j> D (t f ) =o 




,M'/) = arccos 




^(<o) = 0,^(t/) = 0 



(26) 



^D) B U.fc('/) a0 

= 0, = 0 
<M'o) = 0, = o 

Bezeichnungen: 

0 DO Startpunkt Drehwerkswinkel 

0 Df Endpunkt Drehwerkswinkel 

Tla,o Startpunkt Lastposition 

r M ,f Endpunkt Lastposition 

Die Randbedingungen fur den Druck im Zylinder ergeben sich aus den stationaren 
Werten im Start- und Zielpunkt nach Gleichung (11). 

Im Handhebelbetrieb muft dagegen in den Randbedingungen berUcksichtigt wer- 
den, dass die Bewegung nicht aus einer Ruhelage startet und im allgemeinen auch 
nicht in einer Ruhelage endet. Die Randbedingungen zum Startzeitpunkt des Opti- 
mierungshorizonts t Q ergeben sich aus dem aktuellen Systemzustand x(t 0 ), der 
gemessen wird bzw. (iber ein mitgefQhrtes Modell aus den AnsteuergrSIien u SitD , 
u stA und den Medgrofcen <$D>V'D><VA>VA>Pz y l uber einen parameteradaptiven Zu- 
standsbeobachter rekonstruiert wird. 

Die Randbedingungen am Ende des Optimierungshorizonts t f sind frei. 



Aufgrund der technischen Parameter des Kransystems ergeben sich eine Reihe 
von Restriktionen, die unabhangig vom Betriebsmodus im Optimalsteuerungs- 
problem zu beriicksichtigen sind. So ist die Antriebsleistung limitiert. Dies kann Qber 
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einen maximalen Forderstrom in den hydraulischen Antrieben beschrieben werden 
und Gber Amplitudenbeschrankungen fur die Steuergroden in das Optimalsteue- 
rungsprobiem einbezogen werden. 



" U StD,max ~ U StD (0 ~ U StD,max ^ 
' U StA,max ~ U StA ,max 



Zur Vermeidung von Beanspruchungen des Systems durch abrupte Lastwechsel, 
deren Folgen in dem oben beschriebenen vereinfachten dynamischen Modell nicht 
erfasst sind, wird die Anderungsgeschwindigkeit der SteuergroIJen limitiert. Da- 
durch kann definiert die mechanische Beanspruchung limitiert werden. 



* U St D,max — ^StD (0 - ^StD.taax (28) 
m UstA>mx — U StA ,max 



Zusatzlich kann gefordert werden, dass die SteuergroBen als Funktionen der Zeit 
stetig sein sollen und stetige erste Ableitungen bezQglich der Zeit besitzen. 

Der Aufrichtwinke) ist aufgrund der Krankonstruktion limitiert 

<Pa^^<Pa(!)<<P a ^ (29) 
Bezeichnungen: 

UstD,max Maximalwert Ansteuerfunktion Drehwerk 

^stD max Maximale Anderungsgeschwindigkeit Ansteuerfunktion Drehwerk 

ustA,max Maximalwert Ansteuerfunktion Wippwerk 

ti c A Maximale Anderungsgeschwindigkeit Ansteuerfunktion Wippwerk 

St A ,max *"* 

<pA,min Minimalwert Aufrichtwinkel 

<pA,max Maximalwert Aufrichtwinkel 



WO 2004/106215 



PCT/EP2004/005734 



-22- 



Zusatzliche Restriktionen ergeben sich aus weitergehenden Anforderungen an die 
Bewegung der Last. So kann beim vollautomatischen Betrieb, bei dem die gesamte 
Lastbewegung vom Start- bis zum Zielpunkt betrachtet wird, eine monotone Ande- 
rung des Drehwinkels gefordert werden. 



Bahnkorridore lassen sich sowohl im vollautomatischen als auch im Handhebelbe- 
trieb uber die analytische Beschreibung der zulassigen Lastpositionen mit Hilfe von 
Ungleichungsrestriktionen in die Berechnung der optimaien Steuerung einbeziehen. 



Mit Hilfe dieser Ungleichungsbedingungen wird ein Bahnverlauf im Inneren eines 
zulassigen Bereichs, hier des Bahnkorridors erzwungen, die Grenzen dieses zulas- 
sigen Bereichs begrenzen die Lastbewegung und stellen somit ,virtuelle Wande' 
dar. 

Besteht die abzufahrende Bahn nicht nur aus Start- und Zielpunkt, sondern sind 
weitere Punkte in vorgegebener Reihenfolge abzufahren, so kann das durch innere 
Randbedingungen in das Optimalsteuerungsproblem einbezogen werden. 



<P D (t)(<p D (t f )-<p D (t o ))>0 



(30) 



& min 



max 



(31) 



\ 



<Pd (fi ) = <Pd, i > <Pa Ci ) = arocos 



(32) 



Bezeichnungen: 
ft 



(freier) Zeitpunkt des Erreichens des vorgegebenen Bahnpunktes i 
Drehwinkelkoordinate des vorgegebenen Bahnpunktes i 
Radiale Position des vorgegebenen Bahnpunktes i 
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Der Anspruch ist nicht an eine bestimmte Methode zur numerischen Berechnung 
der optimalen Steuerungen gebunden. Der Anspruch bezieht sich ausdrucklich 
auch auf eine angenaherte Losung der oben angegebenen Optimalsteuerungs- 
probleme, bei der im Hinblick auf einen reduzierten Rechenaufwand beim on-line 
Einsatz nur eine Losung ausreichender (nicht maximaler) Genauigkeit ermittelt 
wird. Zudem kann aus Effektivitatsgrunden eine Reihe der oben formulierten harten 
Beschrankungen (Randbedingungen Oder Trajektorienungleichungsbeschrankun- 
gen) numerisch als weiche Beschrankung tiber eine Bewertung der BeschrSn- 
kungsverletzung im Zielfunktional behandelt werden. 

Beispielhaft soil hier jedoch die numerische Losung mittels Mehrstufen- 
Steuerungsparametrisierung eriSutert werden. 

Zur nSherungsweisen numerischen Losung des Optimalsteuerungsproblems wird 
der Optimierungshorizont diskretisiert. 

t 0 =t° <t* <~-<t K = t f (33): 

Die Lange der Teilintervalle [f.f* 1 ] kann dabei an die Dynamik des Problems ange- 
passt werden. Eine grofcere Anzahl von Teilintervallen fuhrt in der Regel zu einer 
Verbesserung der NaherungslOsung, aber auch zu einem erhohten Berechnungs- 
aufwand. 

Auf jedem dieser Teilintervalle wird nun der Zeitverlauf der SteuergrSBen durch ei- 
ne Ansatzfunktion if mit einer festen Anzahl von Parametern (Steuerungspara- 
meter) approximiert. 

u(t)~u apP (0 = U k {t,u k ), t k £t£t M (34) 

Nun kann die Zustandsdifferentialgleichung des dynamischen Modells numerisch 
integriert und das Zielfunktional ausgewertet werden, wobei anstelle der Steuergro- 
Ben die approximierten ZeitverlSufe eingesetzt werden. Im Ergebnis wird das Ziel- 
funktional als Funktion der Steuerungsparameter if, k=0,...,K-1 erhalten. Die 
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Randbedingungen und die TrajektorienbeschrSnkungen lassen sich ebenso als 
Funktionen der Steuerungsparameterauffassen. 

Das Optimalsteuerungsproblem wird auf diese Weise durch ein nichtlineares Opti- 
mierungsproblem in den Steuerungsparametern angenahert, wobei Zielfunktionsbe- 
rechnung und Beschrankungsauswertung des nichtlinearen Optimierungsproblems 
jeweils die numerische Integration des dynamischen Modelis unter BerQcksichti- 
gung des Approximationsansatzes nach Gleichung (34) erfordern. 
Dieses beschrankte nichtlineare Optimierungsproblem kann nun numerisch gelost 
werden, wobei hierzu ein ubliches Verfahren der sequentiellen quadratischen Pro- 
grammierung (SQP) eingesetzt wird, bei dem die LGsung des nichtlinearen Prob- 
lems Ober eine Folge von linear-quadratischen NSherungen bestimmt wird. 
Die Effizienz der numerischen Losung kann betrachtlich gesteigert werden, wenn 
zusatzlich zu den Steuerungsparametern des Intervalls k auch noch der Anfangs- 
zustand 

x k **x(t k \ k = 0,...,K (35) 

des jeweiligen Intervalls als Variable des nichtlinearen Optimierungsproblems be- 
trachtet wird. Durch geeignete GleichungsbeschrSnkungen ist die Stetigkeit der 
approximierten Zustandstrajektorien zu sichem. Damit steigt die Dimension des 
nichtlinearen Optimierungsproblems an. Es ergibt sich jedoch eine betrachtliche 
Vereinfachung in der Verkopplung der Problemvariablen und zudem eine starke 
Strukturierung des nichtlinearen Optimierungsproblems. Daher sinkt der Losungs- 
aufwand in vielen Fallen betrachtlich, vorausgesetzt, die Problemstruktur wird im 
Losungsalgorithmus geeignet ausgenutzt. 

Eine zusatzliche deutliche Reduktion des Rechenaufwands zur Losung des Opti- 
malsteuerungsproblems wird durch eine Approximation mittels Linearisierung der 
Systemgleichungen erreicht. Dabei werden die ursprunglichen nichtlinearen Zu- 
standsdifferentialgleichungen und algebraischen Ausgangsgleichungen (20) ent- 
lang einer zunachst beliebig vorgegebenen Systemtrajektorie (x mf (t), u mf (t)), die die 
Zustandsdifferentialgleichungen erRllIt, linearisiert. 
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Ay = C(t)Ax 

Dabei bezeichnen die GrblJen Ax, Aw, Ay die Abweichungen vom Referenzveriauf 
der jeweiligen GrOBe 

Ax = x-x re/ , Au = M-w re/ , Ay = y-y re/ 

Kef = a ( X ref) + • U re/ (37) 

Die zeitvarianten Matrizen Aft), B(t), C(t) ergeben sich aus den Jacobi-Matrizen 



Wird nun die Optimalsteuerungsaufgabe in den Variablen 

Ax,A« formuliert, so ergibt sich ein beschranktes linear-quadratisches Optimal- 
steuerungsproblem. Bei geeigneter Wahl der Ansatzfunktionen if kann die Zu- 
standsd'rfferentialgleichung uber die zugehorige Bewegungsgleichung auf jedem 
Teilintervall if,?* 1 ] analytisch gelbst werden, und die aufwendige numerische Integ- 
ration entfallt. 

Die Optimalsteuerungsaufgabe wird so durch ein endlich-dimensionales quadrati- 
sches Optimierungsproblem mit linearen Gleichungs- und Ungleichungsrestriktio- 
nen approximiert, welches mit einem angepassten Standardverfahren numerisch 
gelbst werden kann. Der numerische Aufwand hierftir ist wiederum deutlich gerin- 
ger als bei dem oben beschriebenen nichtlinearen Optimierungsproblem. 
Der beschriebene Linearisierungsansatz ist besonders fur die naherungsweise L6- 
sung der Optimalsteuerungsprobleme bei Handhebelbetrieb geeignet, da in diesem 
Fall zum einen aufgrund des kOrzeren Optimierungshorizonts (Zeitfenster [t 0 J f ]) 
die durch die Linearisierung bedingten Ungenauigkeiten eine geringere Auswirkung 
haben und zum anderen mit den im jeweils vorangehenden Zeitschritt berechneten 



A(t) = 



d(a(x ref {t)) + b{x ref (t))u^{t)) 
dx^if) 



B(t) = b(x re At)), C(0 = 



dc(x nf (t)) 



..(38) 
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optimalen Steuerungs- und Zustandsverlaufen geeignete Referenztrajektorien zur 
Verfiigung stehen. 

Als Losung des Optimalsteuerungsproblems werden die optimalen Zeitverlaufe so- 
wohl der Steuergr6R>en als auch der Zustandsgroften des dynamischen Modells 
erhalten. Diese werden bei Betrieb mit unterlagerter Regeiung als Stell- und Fuh- 
rungsgroden aufgeschaltet. Da in diesen Sollfunktionen das dynamische Verhalten 
des Krans berucksichtigt ist, miissen durch die Regeiung nur noch StorgrSften und 
Modellabweichungen ausgeglichen werden. 

Bei Betrieb ohne unterlagerte Regeiung werden die optimalen Veriaufe der Steuer- 
grofien dagegen direkt als StellgroBen aufgeschaltet. 

Weiterhin liefert die Losung des Optimalsteuerungsproblems eine Prognose der 
Bahn der pendelnden Last, die fur erweiterte Madnahmen zur Kollisionsvermeidung 
nutzbarist. 

Fig. 8 zeigt das Ablaufdiagramm fGr die Berechnung der optimierten Steuergrofce 
im vollautomatischen Betrieb. Dies untersetzt Modul 37 aus Fig. 3. Ausgehend von 
den durch die Sollpunktmatrix festgelegten Start- und Zielpunkten der Lastbewe- 
gung wird das Optimalsteuerungsproblem durch Einbeziehung der Vorgabe des 
zulassigen Bereichs und der technischen Parameter definiert. Die numerische Lo- 
sung des Optimalsteuerungsproblems liefert optimale Zeitverlaufe der Steuer- und 
Zustandsgroden. Diese werden bei unterlagerter Regeiung zur Lastpendeldamp- 
fung als Stell- und Fuhrungsgrolien aufgeschaltet. Alternativ kann eine Realisierung 
ohne unterlagerte Regeiung - dann mit direkter Aufschaltung der optimalen An- 
steuerfunktionen auf die Hydraulik - realisiert werden. 

Fig. 9 zeigt das Zusammenwirken von Zustandsrekonstruktion und Berechnung der 
optimalen Steuerung im Fall des Handhebelbetriebs. Der Zustand des dynami- 
schen Kranmodells wird unter Nutzung der verfOgbaren MeligroRen nachgefQhrt. 
Durch Losung des Optimalsteuerungsproblems werden solche Zeitverlaufe der An- 
steuerfunktionen ermittelt, die - ausgehend von diesem akuellen Zustand - bei re- 
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duziertem Lastpendeln die Lastgeschwindigkeit an die uber die Handhebel vorge- 
gebenen Sollwerte heranfdhrt. 

Eine einmal berechnete optimale Steuerung wird nicht uber den vollen Zeithorizont 
[t 0 J f ] realisiert, sondern laufend an den aktuellen Systemzustand und die aktuel- 

len Sollwerte angepasst. Die Frequenz dieser Anpassung ist durch die benotigte 
Rechenzeit zur Neuberechnung der optimalen Steuerung limitiert. 

Fig. 10 zeigt beispielhafte Ergebnisse fur optimale Zeitverlaufe der Steuergrdften im 
vollautomatischen Betrieb. Dabei wurde ein Zeithorizont von 30 s vorgegeben. Die 
Ansteuerfunktionen sind stetige Funktionen der Zeit mit stetigen 1. Ableitungen. 

Fig. 1 1 zeigt beispielhafte Zeitverlaufe von Steuergrofcen und RegelgrofJen im si- 
muiierten Handhebelbetrieb. Die Sollwerte fur die Lastgeschwindigkeit (die Hand- 
hebelvorgaben) werden in Form von zeitlich versetzten Rechteckimpulsen variiert. 
Die Aktualisierung der optimalen Steuerung erfolgt mit einer Abtastzeit von 0.2 s. 
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Kran Oder Bagger zum Umschlagen von einer an einem Lastseil 
hangenden Last mit optimierter Bewegungsfiihrung 



Patentanspriiche 



1 . Kran oder Bagger zum Umschlagen von einer an einem Lastseil hangenden 
Last mit einem Drehwerk zum Drehen des Kranes oder Baggers, einem 
Wippwerk zum Aufrichten bzw. Neigen eines Auslegers und einem Hubwerk 
zum Heben bzw. Senken der an dem Seil aufgehangten Last mit einem An- 
triebssystem, 

gekennzeichnet durch 

eine Bahnsteuerung (31) deren AusgangsgrofJen (u 0U tD, u ou tA. u ou tL, u 0 uir) direkt 
oder indirekt ais Eingangsgr6R>en in die Regelung fur die Position bzw. die 
Geschwindigkeit des Kranes (41) oder Baggers eingehen, wobei die Ftih- 
rungsgrSften for die Steuerung in der Bahnsteuerung (31) so generiert wer- 
den, daB sich eine Lastbewegung mit minimierten Pendelausschlagen ergibt. 
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2. Kran oder Bagger nach Anspruch 1, dadurch gekennzeichnet dad in der 
Bahnsteuerung (31) modellbasiert optimale Steuertrajektorien online berech- 
net und aktualisiert werden. 

3. Kran oder Bagger nach Anspruch 2 gekennzeichnet durch modellbasierte op- 
timale Steuertrajektorien basierend auf einem um Referenztrajektorien lineari- 
siertem Modell. 

4. Kran oder Bagger nach Anspruch 2 gekennzeichnet durch modellbasierte op- 
timale Steuertrajektorien basierend auf einem nichtlinearen Modellansatz. 

5. Kran oder Bagger nach einem der Anspruche 1 bis 4 gekennzeichnet durch 
modellbasierte optimale Steuertrajektorien unter RQckfQhrung von alien Zu- 
standsgroBen. 

6. Kran oder Bagger nach einem der Anspruche 1 bis 4 gekennzeichnet durch 
modellbasierte optimale Steuertrajektorien unter Ruckfuhrung von mindestens 
einer MeRrofte und Schatzung der verbleibenden ZustandsgroRen. 

7. Kran oder Bagger nach einem der Anspruche 1 bis 4 gekennzeichnet durch 
modellbasierte optimale Steuertrajektorien unter RuckfQhrung von mindestens 
einer Mefcrofce und Nachfuhrung der verbleibenden ZustandsgroRen durch 
modellbasierte Vorwartssteuerung. 

8. Kran oder Bagger nach einem der Anspruche 1 bis 7, dadurch gekennzeich- 
net, daft die Bahnsteuerung (31) als Vollautomatik oder als Halbautomatik 
durchfQhrbar ist. 

9. Kran oder Bagger nach einem der vorangehenden Anspruche, dadurch ge- 
kennzeichnet, dad in die Bahnsteuerung (31) als Eingangsgrofte eine Soll- 
punktmatrix (35) fur die Position und Orientierung der Last eingebbar ist. 
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10. Kran oder Bagger nach einem der vorangehenden Anspriiche, dadurch ge- 
kennzeichnet, daft die Sollpunktmatrix (35) aus Start- und Zielpunkt besteht. 

11. Kran oder Bagger einem der vorangehenden Anspriiche, dadurch gekenn- 
zeichnet, daft im Fall eines halbautomatischen Betriebs zusatzlich die ge- 
wunschte Zielgeschwindigkeit der Last durch die Stellung des Handhebels 
(34) in die Bahnsteuerung (31 ) eingebbar ist. 

12. Kran oder Bagger nach Anspruch 11, dadurch gekennzeichnet, daft im halb- 
automatischen Betrieb uber Sensoren die Meftgroften der Positionen von 
Kran und Last erfaftbar und in die Bahnsteuerung (31 ) zuriickfuhrbar sind. 

13. Kran oder Bagger nach Anspruch 11, dadurch gekennzeichnet, daft im halb- 
automatischen Betrieb die Positionen von Kran und Last in einem Modul zum 
modellbasierten Schatzverfahren (43) schatzbar und in die Bahnsteuerung 
(31 ) zuruckfQhrbar sind. 

14. Kran oder Bagger nach einem der vorangehenden Anspriiche, dadurch 
gekennzeichnet, daft die Ausgangsgroften (u 0 utD, u ou tA, u ou tL. u out R) zunachst 
in eine unterlagerte Regelung mit Lastpendeldamfung gefuhrt werden. 

15. Kran oder Bagger nach Anspruch 14, dadurch gekennzeichnet, daft die 
LastpendeldSmpfung ein Bahnplanungsmodul, eine Zentripetalkraftkompen- 
sationseinrichtung und zumindest einen Achsregler fur das Drehwerk, einen 
Achsregler fur das Wippwerk, einen Achsregler fur das Hubwerk und einen 
Achsregler fQr das Schwenkwerk aufweist 

16. Kran oder Bagger nach einem der vorangehenden Anspriiche, dadurch 
gekennzeichnet, daft mittels der Bahnsteuerung (31) die Bewegungs- 
bahn der Last derart festlegbar ist, daft vorbestimmte freie Bereiche von 
der pendelnden Last nicht verlassen werden konnen. 
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Fig.2 
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Fig.3 
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Fig.4 
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Fig. 8 
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Fig. 10 
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